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ABSTRACT 



Aims. We investigate the reconstruction of pulsar magnetospheres close to the light cylinder surface to study the curvature drift 
instability (GDI) responsible for the twisting of magnetic field lines in the mentioned zone. The influence of plasma dynamics on the 
saturation process of the GDI is studied. 

Methods. On the basis of the Euler, continuity, and induction equations, we derive the increment of the GDI and analyze parametric ally 
excited drift modes. The dynamics of the reconstruction of the pulsar magnetosphere is studied analytically. 

Results. We show that there is a possibility of a parametrically excited rotational-energy pumping-process in the drift modes. It is 
indicated by the generation of a toroidal component of the magnetic field that transforms the field lines into such a configuration, in 
which plasma particles do not experience any forces. At this stage, the instability process saturates and the further amplification of 
the toroidal component to the magnetic field lines is suspended. 
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1. Introduction 

Exploring the magnetic field of the Crab nebula, 
[Piddington (1957)] was the first to discover the presence of 
a central object in the nebula, with frozen magnetic field inside. 
It was assumed that the rotation of a central object provokes 
the generation of the toroidal component of magnetic field. 
Investigations have also shown that this kind of magnetic field 
characterizes magnetized star outflows ( IWeber & Davis 1967I ). 
Pulsars are one of the most distinctive examples of rotationally 
powered magnetized stars producing the relativistic outflows 
commonly known as pulsar winds. 

One of the major problems for pulsar winds concerns the 
transition of the magnetized plasma flows by means of the so- 
called "light cylinder" surface, a hypothetical zone, where the 
linear velocity of rotation equals the speed of light. According to 
the model of pulsar magnetosphere jGoldreich & Julian 19691 ). 
the relativistic plasma flow, emanating from the pulsar sur- 
face streams along very strong magnetic field lines. The typ- 
ical values of the magnetic field are of the order of 10'^ 
G (Manchester (& Taylor 19771. Such a strong magnetic field 



forces the particles, streaming toward the light cylinder zone, 
to co-rotate with a pulsar On the other hand, once the rela- 
tivistic effect of the mass increment is taken into account, the 
radial acceleration of the particles appears to be limited. This 
means that the particles will never cross the light cylinder surface 



if the rigid rotation is preserved ( Machabeli & Rogava 19941. 



Therefore, the co-rotation of the magnetized plasma cannot be 
maintained close to the mentioned zone and a further exploration 
of this phenomenon is needed. 
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The simplest means of explaining the transition in the plasma 
outflow at the light cylinder surface is to constrain the mo- 
tion of the flow to become asymptotically close to a regime, 
in which particles do not experience any forces (see e.g., 
Machab eli et al. 20001 [Rogava et al. 2003| l. For this purpose, if 
the magnetic field is still strong, one needs to twist the field 
lines in an appropriate way. In other words, it is necessary to 
generate the toroidal component of the magnetic field. The pul- 
sars maintain the dipolar magnetic field with a certain curvature 
of the field lines. The particles moving along the curved field 
lines drift perpendicularly to the curvature plane. The existence 
of the drift motion creates the necessary conditions for develop 
ing a GDI (see e.g., [Kazbegi et al. 1989| [Kazbegi et al. 1991a 



Kazbegi et al. 1991b Shapakidze et. al 2003) l. The mechanism 
excites transverse drift modes providing the magnetic field with 
the toroidal component. The GDI becomes effective in the vicin- 
ity of the light cylinder surface, where the pulsar's dipolar field 
is comparable to the toroidal component of the excited mode. 
The kinetic energy of the primary beam particles, of Goldreich- 
Julian (GJ) density, can indeed supply the generation of the drift 
waves. Mostly, it is assumed that the toroidal component of the 
magnetic field is created by the GJ current. However, the kinetic 
energy density, Wb, of the primary beam is insufficient to change 
the configuration of the dipolar magnetic field, Bq, significantly, 
since Wb « B^/{47:). Therefore, the process of generating the 
toroidal component, B,- (of the order of Bo), requires an addi- 
tional energy supply, which can be extracted from the pulsar ro- 
tational energy. We have akeady shown that the rotational en- 
ergy pumping into the drift modes can be implemented by the 
"parametric" instability (lOsmanov et al. 20081 1. This mechanism 
is called parametric, because the effect is caused by the relativis- 
tic centrifugal force, which, as a parameter, changes in time and 
induces the instability. The aim of the present paper is to study 
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Fig. 1. a) The geometry in which the set of main Eqs. ( IH3b is consid- 
ered; (tg, Cr, 6;,-) denotes orthonormal basis of unit vectors; e, is directed 
perpendicularly to the plane of the figure; C is the center of the curva- 
ture; and Cfl and e, are the tangential and perpendicular (with respect 
to the field line) unit vectors respectively, b) The geometry for deriving 
Eq. J23t : the curved line denotes the twisted magnetic field, B, gener- 
ated due to the raising of magnetic perturbation, B, . Note, that B, and 
Be are oriented with respect to the initial 'quasi straight' magnetic field 
line. 



the saturation process of this instability close to the light cylinder 
zone. 

The problem analysis is completed well if one considers the 
mechanical analog of the motion of plasma flow along the pulsar 
magnetic field lines and studies the dynamics of a single particle 
moving inside the rotating channels. In the present paper, we ap- 
ply the method developed in ( Rogava et al. 2003) ) and study the 
plasma process that converts the configuration of the magnetic 
field lines to the Archimedes' spiral: it appears that in order to 
suppress the reaction force, the field lines, rather being straight, 
should deviate back and lag behind the rotation; consequently, 
the curvature will eventually increase, inevitably leading to the 
decrement of the reaction force; this process will last until the 
magnetic field lines form the shape of the Archimedes' spiral, 
which removes the reaction force and saturates the instability. 

The paper is arranged as follows. In section 2, we examine 
the parametric curvature drift instability and derive an expres- 
sion for its corresponding growth rate. In section 3, we consider 
the kinematics of a particle moving along the rotating magnetic 
field line. In section 4, we apply the method for 1 -second pulsars 
and the Crab pulsar, and in section 5 we summarize our results. 



2. Parametric curvature drift instability 

It is well known that the presence of an external varying pa- 
rameter usually generates the plasma instability. For example, 
the mechanism of energy pumping process from the external 
alternating electric field into the electron-ion plasma is quite 
well investigated (see e.g.,|SiUn 19731 paleev & Sagdeev 197J 



IMax I973l l. Although the physics of the parametric instability in 
the electron-ion plasma differs from that of the electron-positron 
(e~e^) plasma, the techniques of calculation can be the same. 
The case of a e^e^ plasma, where the external varying param- 
eter is the altering centrifugal acceleration, was considered by 
Machabeli et al. (2005). 

We start our consideration by initially assuming that the 
magnetic field lines are almost straight with very small non 
zero curvature (see Fig. [T^). The plasma stream is supposed to 
move along the co-rotating field lines of a very strong mag- 
netic field (~ 10'^ G for typical pulsars). We assume that the 
plasma flow consists of two components: the plasma component 
composed of electrons and positrons (e^); and, the so-called, 
beam component (b) composed of ultra relativistic electrons. 
The dynamics of plasma particles moving along the straight ro- 



tating magnetic field lines is described by the Euler equation 
(IMachabeli et al. 20051 1 



+ (v„V)p„ 

ot 



where 



-c^r«^V^ +-|E + -v„xB 
m \ c 



(1) 



R is the coordinate along the straight field lines, a = {e'^,b} 
denotes the sort of particles, p^,, Vq,, and ea are the momentum 
(normalized to the particle's mass), the velocity, and the charge 
of electrons/positrons, respectively; E is the electric field and B 
is the magnetic field. On the right-hand side of Eq. ([1]), there 
are two major terms, the first of which represents the centrifugal 
force and the second, the Lorentz force. The full set of equations 
for n, V, E, and B should be completed by the continuity equation 



— + "^(riaVa) = 0, 
Ot 

and the induction equation 
1 5E 47r 

z 



V xB = -— + 

c ot c 



Zj^' 



(2) 



(3) 



where «„ and Jq. are the density and the current, respectively. 

Rewriting the Euler equation in Eq. ([T) for the leading state 
and taking into account the frozen-in condition Eq + ^\oo x Bq = 
0, the solution for ultra relativistic particle velocities writes as 
follows ( [Machabeli & Rogava 1994| i 



: C COS(Qf + if), 



(4) 



where Vj denotes the velocity component along the magnetic 
field lines and if, the initial phase of each particle. 

The centrifugal force eventually causes the separation of 
charges in plasma consisting of several species. This process be- 
comes so important that the coiTesponding electromagnetic field 
affects the dynamics of the charged particles. Therefore, the pro- 
duced electric field should also be considered in Eq. ([TJ as the 
next term to be approximated (Osmanov et al. 2008). 

To clarify the mathematical treatment, we linearize the set 
of Eqs. ([T][3]l, assuming that, to the zeroth approximation, the 
flow has a longitudinal velocity satisfying Eq. (HJl and also drifts 
along the x-axis because of the curvature of magnetic field lines 
(see Fig.[T]i: 



Uq, — 



<^b„Rb 



(5) 



where m„ is the drift velocity; ojb^ - eaBo/mc, Rb is the cur- 
vature radius of magnetic field lines, and Bo is the magnetic in- 
duction in the leading state. In our case, the curvature drift con- 
tributes to J [see Eq. Q] as a source of additional cuiTent. 

We represent all physical quantities as the sum of the zeroth 
and the first order terms 



where 



(6) 



= {n,\,p,E,B}. 
We then express the perturbed quantities as 
»I''(f,r)oc>I'i(f)exp[/(kr)] , 



(7) 
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examining only the x components of Eqs. ( |l|3t , considering the 
perturbations with kg <k and kr = 0, and bearing in mind that 
a: cE\/B{), one can easily reduce the set of Eqs. ([T][3l) to the 
form 

dt 



Z E 



Js(ga)Jn(h)Jt(ga)Jp(h) 



II 1-1-7/^ II 



- i(k^Ua + kgv^^nl = ik^nlvl , 
- ikgcBl = 47r ^ ^^.(n" v^^ + n^M„). 

Introducing a special ansatz for and n]^ 
where 

Aa^it) = ^ (Qf + ^) + ^ sin[2 (fif + 0], 
A„„(f) = ^ sin(Qf + ^p), 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
(14) 



ch^,b ^y^o^ec ^jrii (cj+!^+ n(2[s +ju]+ n)) 

X Br (W + Q (2[i -l+fl]+n- p+0-))X e'V(2[.v-/+//]+«-p+<r)^ ^jg-, 

where 



h - — 



To solve Eq. ( fT9] l, we must examine similar equations, e.g., 
rewrite Eq. ( fT9] l for Z?,(w ± O), Z?r(w + 2Q), etc.. This means that 
one has to solve a system with an infinite number of equations, 
which makes the problem impossible to handle. Therefore, the 
only solution is to consider the physics close to the resonance 
condition, which provides the cutoff to the infinite row in Eq. 
( fT9] l making the problem solvable ( ISilin & Tikhonchuk 19701 ). 

Let us consider the resonance condition, which corresponds 
to the curvature drift modes. As is clear from Eq. (fT9T l. the proper 
frequency for the CDI equals 



(20) 



and substituting Eqs. (II 1112b into Eqs. (I8I9I I. one derives the ex- 
pressions 



1 _ t^a ;kA„(0 
mcjao 



j' e-''^^(''\it')Br(t')dt', (15) 



Therefore, physically meaningful solutions relate the case, when 
ki-Ua/2 < 0. The present condition implies that 2i + n = and 
2[s + fi] + n - 0. On the other hand, it is easy to check that, for 
the typical quantities of 1 -ieconc/ pulsars, y^, ~ 10^ and A ~ 6 x 
10'" - 3 X 10"cm (where A ^ - l/k^ is the wavelength), one 
has \k^Ua/2\ ~ (0.02h-0.1).s"' (here, it is assumed that kj, < and 
Uh > 0, otherwise the resonance frequency would be negative). 

Since particles have different phases, to solve Eq. ( fT9] l, we 
must examine the average value of Br with respect to (p. Then, 
by taking into account the formula 



-'^•'^•"\.it")Br{t")dt" . 



(16) 



Combining Eqs. ( I15I16I I with Eq. ( fTOl l. it is straightforward to 
reduce it to the form 



preserving only the leading terms of Eq. ( |T9] l, and also taking 
into account that the beam components exceed the correspond- 
ing plasma terms by many orders of magnitude, one derives the 
dispersion relation for the CDI (lOsmanov et al. 20081 1 



-ikgcBlit)^ y i:ie'"^"« fe-'^'^'\(t')Brit')dt'+ L + 
y i^fe^M^e'"^"*" f'dt' [' e-''^"('">v„(f")B.(f")*",(17) 



2ybokgc 



Jo 



kxUh\ J (kgc 
4Q r" Tl" 



(21) 



To determine the CDI growth rate, F, let us write o) = ll)o + iF 
and substitute this into Eq. dSTT i. Then, it is easy to show that the 
increment is given by 



Where Ua - e s/A-mi^/m is the plasma frequency. To simplify 
Eq. ( fTTI l, one may use the following identity 



(18) 



3 (^l k^ut 

2 Jha kgc ^ 



1/2 



kxUb\ . (kgc 



(22) 



where J six) (s = 0; +1; +2 . . .) is the Bessel function of integer 
order ( Abramovitz & Stegan 1965 1. Eq. ([TtT i then reduces to the 
form 



B 



J s{ga)Jnih)J l{ga)J i,(h) ^ 



xBr {(jj + Q. {2[s -l]+n-p + a-)) 



kxU(y 



co+^+ n(2s + n) 



We can qualitatively analyze how the shape of magnetic field 
lines changes with time. After perturbing the magnetic field in 
the transverse direction, the toroidal component will strengthen 
and the field lines will gradually lag behind the rotation. This 
process will inevitably influence the dynamics of the particles. 
As a result, the acceleration will be lower than in the case of 
straight field lines. Since the CDI is centrifugally excited, the 
corresponding efficiency will also decrease, completely vanish- 
ing for such a configuration of the magnetic field lines, when the 
particles do not accelerate centrifugally at all. 

Let us assume that the corresponding critical value of the 
toroidal component, when the instability diminishes, is B, . Then, 
from simple considerations, we can estimate the corresponding 



4 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 




R/R 



Fig. 2. The arm of Archimedes spiral in polar coordinates (O, R). 
Two orthonormal bases are considered: i) polar components of 
unit vectors, (eo,ejj); and ii) normal and tangential components 
of unit vectors, (e,-, ee), respectively. C is the center of the curva- 
ture. 



timescale of this transformation. Indeed, referring to Fig. [T]3, 
tan;^' = Br/Bg. On the other hand, the toroidal component of 
magnetic field behaves with time as 

Br « By, 

where B° denotes the initial perturbed value of the toroidal com- 
ponent. Consequently, the timescale of the process can be esti- 
mated to be 



1 



In 



r(^0 \Bgtanx 



(23) 



3. Kinematics of particles moving along the curved 
magnetic field lines 

We describe the kinematics of particles moving along the coro- 
tating curved magnetic field lines. In particular, we adopt the 
Archimedes' spiral for the configuration of magnetic field lines 



0) = aR, 



where O and R are the polar coordinates and a - const (see 
Fig. |2]i. Rogava et al. (2003) showed that the dynamics of the 
particle motion asymptotically tends to the force-free regime, if 
the particle slides along the rotating channel that has the shape of 
the Archimedes' spiral. However, in the laboratory frame (LF), 
the particle follows straight linear paths with constant velocities. 
In this case, an observer from the LF will measure the effective 
angular velocity 



Q-,f = Q H — D. + flv. 



(25) 



where v is the radial velocity of a particle motion, and Q. is the 
angular velocity of rotation. If the particle moves without accel- 
eration in the LF, the corresponding effective angular velocity 



Fig. 3. Behaviour of /3r versus R/Rc- The set of parameters is 
p = Is,/3rq = {0.01; 0.3; 0.5; 0.7; 0.9). 



must be equal to zero (Of// = 0). In this case, Eq. 
that 



O 

Vc = , 

a 



infers 



(26) 



which means that, for any Archimedes' spiral with a < -Q./c, 
the LF trajectory of the particle becomes a straight line if it 
moves with a certain "characteristic velocity", Vc- 
The relativistic momentum of the particle is given by 



Pr — ymv, 
- ymRQef, 



(27) 
(28) 



where m and y are the rest mass and the Lorentz factor of the 
particle, respectively. The expansion of the equation of motion, 
dP/dt - F {F is the reaction force), in terms of the radial com- 
ponent 



aR 



:|F|, 



(24) = 



and the angular component 
1 



(29) 



(30) 



of the reaction force, yields the equation (Rogava et al. 2003 1 

d^R Q. - y^v{a + ^Iv/c^) 



dfi 

where 



-Q.efR, 



Q^R^ 



(31) 



a^R^ 



1/2 



We note that {d^Rjdf-) = 0, when v = Vc = -Q/a. 

Since we are interested in relativistic flows, let us consider 
Vc - c setting a - -Q./c. In Fig. |3] we plot the solution to Eq. 
(ISTT i. that is the dependence of /3r = v/c versus R/Rc (Rc is the 
light cylinder radius) for the different initial values of /3ro. As 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



Fig. 4. a) The trajectory of a particle in the Rotational Frame (RF) of 
reference; and b) the trajectory of a particle in the Laboratory Frame 
(LF) of reference. The set of parameters is P = Is and jSj^^ = 0.01. The 
radial distances are taken in terms of the light cylinder radius, R/Rc- 
The trajectory of the particle in RF follows the Archimedes' spiral (a), 
while the trajectory in the LF asymptotically tends to a straight line 
configuration (b). 



is clear from these plots, even though the initial velocity of par- 
ticles is weakly relativistic (e.g., (Sro = 0.01), it asymptotically 
converges to the characteristic velocity, v^. As a result, the parti- 
cle trajectory in the LF must become linear. 

Indeed, in Fig. |4] we show the particle trajectory in the 
Rotational Frame (RF) of reference (see Fig. |4^) as well as in 
the LF (see Fig. |4}3) for both /3ro - 0.01 and the same spiral 
configuration with a = -Q/c. Observing the particle trajectory 
from the LF, one can note that the path asymptotically becomes 
linear, indicating that the rotational energy pumping process di- 
minishes. Therefore, we conclude that the magnetic field with 
the of the Archimedes' spiral may guarantee the saturation pro- 
cess of the CDL 

On the other hand, the equation of Archimedes' spiral yields 
tan;^' - aR - Br/Bg (see Fig.[T}5), resulting in B, = Bg at the LC 
when a = -Q/c. 

4. Discussion 

We consider Eq. ( |23T l and plot the timescale of the instability ver- 
sus the wavelength, /} for several values of initial toroidal per- 
turbations B*/. In Fig.|5] we display the behavior of T{A^), when 
tanx = 1, for several values of initial perturbation B^J Bg = 5 e 
{10"'; lO"-'; lO"''; 10"^). Two major applications are examined; 
(a) 1 -second pulsar and (b) the Crab pulsar From Eq. ( |23] ). we 
can infer that the timescale is a continuously decreasing function 
of the initial perturbation, Bf; as becomes smaller the S param- 
eter, so does the initial perturbation and, consequently, the mag- 
netic field lines need more time to achieve the required struc- 
ture. As observations show, the ratio PIP ranges from 10" 
(PSR 0531 - Crab pulsar) to 10'^? (PSR 1952+29). However, 
the greatest values of the twisting timescales, shown in Fig. |5] 
vary between lO^'i and 10^.?, which are shorter by many orders 
of magnitude than Pj P, illustrating the high efficiency of the in- 
stability. 

The reconstruction of the magnetic field requires a certain 
amount of energy, therefore it is essential to estimate pulsar's 
slowdown luminosity (Lp) and compare it to the, so-called, 
"magnetic luminosity" (L,„ = AEm/T, where aE„ is the vari- 
ation in the magnetic field energy due to the instability). For Lp, 
one has 



47T-P 



(32) 



Fig. 5. Here we show the behaviour of the timescale of sweepback of 
the magnetic field lines versus the wave length, A « /!, , for the 1 -second 
pulsar [see (a)] and for the Crab pulsar, P « 0.033.S [see (b)], respec- 
tively. The set of parameters is: ^ = {10"'; 10"^ 10"^; 10"'), yj, ~ 10* 
andifl = lOOOi?,:- 



respectively. For typical pulsars, PjP 
luminosity is estimated to be 

Lp K 7.9 X 10^'erg/i. 



, the slowdown 



(33) 



The "magnetic luminosity", L„„ can be estimated straightfor- 
wardly 

L„, = -^AV, (34) 

where aV ~ R^Ar - RIs (e = Ar/Rc « 1) is the volume, 
where the process of twisting takes place. Since we are study- 
ing the instability in the nearby zone of the LC, the magnetic 
field components must be given as B^) ~ 6 ■ Bg - 6 ■ BpiRplRcf , 
where Bp ~ lO'^G is the magnetic field near the pulsar's surface. 
Substituting all quantities into Eq. ( l34l l corresponding to the en- 
ergy gain for A ^ ~ 2.5 x 10"cot, and T « 10^ s, (see Fig.|5^), and 
bearing in mind Eqs. ( I2I23I ), one can show that for e ~ 0.1 and 
6 ~ Q.l the value of the "magnetic luminosity" is approximately 
equal to 



L,„ ^ 7.3 X lO^^erg/s 



(35) 



lC™^ ^ 4 7 X lO^^erg/s. 



where / ~ MRJ, is the moment of inertia of the pulsar, and M ~ 
Mq X 2x lO-'^'g and Rp ~ lO^cm are the pulsar's mass and radius. 



A direct comparison between Eqs. ( [33] ) and ( [35] l infers that 
L„, << Lp, meaning that approximately only 0.0001% of the 
total energy budget goes to the reconstruction of the magnetic 
field lines. 

In Fig.lSj), we show the behavior in T{A^) for the Crab pulsar 
in the different cases of initial magnetic perturbations. The de- 
pendence does not change qualitatively but the twisting process 
changes quantitatively as the corresponding timescale is now of 
the order of ~ 10:10^ s. The "magnetic luminosity" of the Crab 
pulsar for 6 = O.i, = 7 x IQ^cm, T x 4s (see Fig.|5]D), and 
e = 0.1 approximately equals [see Eq. ( l34] il 

(36) 

On the other hand, a direct calculation of the Crab pulsar {P ^ 
0.033 s) luminosity, estimated by Eq. (l32l i. yields 

L'p"''' ^ 9.3 X lO^'^erg/s. (37) 

Therefore, as in the case of the Crab pulsar, the energy required 
to reconstruct the magnetosphere averages ~ 0.00005% of the 
pulsar's energy and the twisting process becomes feasible in this 
case as well. 

The sweepback mechanism described in this paper, there- 
fore, appears to be extremely efficient for pulsar magneto- 
spheres. The curvature drift instability may lead to the recon- 
struction of the magnetic field lines in such a way that the dy- 
namics of magnetosphere becomes force-free, which in its turn, 
completely disrupts the instability, thus saturating the CDI. 
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5. Summary 

We have analytically examined the non-stationary pattern of the 
pulsar magnetosphere near the LC zone. The GJ current can be 
estimated by Jcj - etihc, giving the value of the corresponding 
magnetic field, B, * AjiJcjR„Ic, where R„ x BiR)/B'{R) = R/3 
(where B' = dB/dR) is the length scale of the spatial inhomo- 
geneity of the magnetic field. If we assume a dipolar configura- 
tion, then, taking the value of the GJ density, ri/, x Q.B l{2nec), 
into account, one can show that the toroidal magnetic field 
equals -^B. Inside the light cylinder {R < R^), this value is 
less than the required one-Z? and reaches its maximum value, 
|B, on the light surface. Therefore, the GJ current cannot sig- 
nificantly change the configuration of the magnetic field. This 
current evidently exceeds the curvature drift current, because 
Jcj/Jcur ~ c/uh » 1. However, the drift current. Jar, is not 
the source of the toroidal component, B^, of the magnetic field, 
but it is a trigger mechanism for generation of the perturbed 
current, Ji = ^("i^], + n^Mi) responsible for the creation of 
Br (see Eq. (fTOl)). The source of the instability of the current 
and the resulting magnetic field is the pulsar's rotational energy. 
We have found that the instability is achieved by the paramet- 
ric mechanism, which effectively pumps the pulsar rotational 
energy directly into the generated mode. The process lasts un- 
til the plasma dynamics reaches the force-free regime of mo- 
tion and the overall magnetosphere relaxes to the steady state 
configuration. |Bucciantini et al. (2006)| considered the dynam- 
ics of rotating pulsar winds by performing the numerical solu- 
tion of relativistic magnetohydrodynamic (RMHD) equations in 
the Schwarzschild metric. The system was allowed to relax to 
a steady state configuration rapidly approaching the force-free 
regime. We suppose that the full set of RMHD equations com- 
prises the terms applicable to CDI. Therefore, we highlight the 
numerical results obtained by |Bucciantini et al. (2006)| as confir- 
mation of our theory. On the other hand, the aim of the present 
work was to explain the saturation of the CDI in terms of the 
generation of currents, which makes the physics of the process 
more transparent. 

The main aspects of the present work can be summarized as 
follows: 

1 . Examining the pulsar magnetospheric relativistic plasma, we 
have studied the role of the parametrically excited CDI in the 
process of sweepback of magnetic field lines and the satura- 
tion process of the instability. The present parametric insta- 
bility is based on the method developed by Silin (1973), but 
differs from this in principal, since we study an alternating 
centrifugal force instead of an alternating electric force. 

2. The linear analysis of the Euler, continuity, and induction 
equations yields the dispersion relation governing the CDI. 

3. Considering the resonance frequencies of the sweepback 
process, an expression of the instability increment has been 
obtained. 

4. On the basis of the expression of the instability growth rate, 
we have derived the formula of the transition timescale of 
quasi-linear configuration of field lines into the Archimedes' 
spiral. The particles' motion is force-free along these mag- 
netic field lines, leading to the saturation of the instability. 

5. The transition timescale has been compared to the wave- 
length for 1 -second pulsars and the Crab pulsar. For both 
cases it was shown that the corresponding timescale is 
shorter than pulsar's spin down rates by many orders of mag- 
nitude illustrating the high efficiency of the discussed pro- 
cess. 
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